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We investigate elliptic flow in heavy ion collisions at intermediate energies. In doing this, we 
implement and use a lattice-Hamiltonian model of the nuclear interaction and we also study the 
effect of in-medium nucleon-nucleon cross sections that follow consistently from the momentum- 
dependence of the nuclear mean field. 



... I. INTRODUCTION 

(N ■ 

■ Heavy ion collisions is an area of research which seeks to study physical systems under extreme conditions of density 
\ and temperature in experiments occuring in terrestrial laboratories. As such, it represents a rich and challenging field 
. well worthy of intellectual pursuit and its study is important for a deep and a more complete understanding of Nature. 
' At high energies, one goal of this program is to form and study a new state of matter which is a tantalizing prediction 

of QCD: the quark-gluon plasma [Q. At lower energies, the experimental and theoretical efforts have focused on the 
need to characterize and quantify the nuclear equation of state E|. This physics also has an important role to play 

04 in the theory of supernovae and that of neutron star properties It is that energy regime that we consider in this 

O^ ' work. 

\ In order to identify novel many-body features without ambiguity it is imperative to provide a realistic model of 
^ ■ the nuclear reaction dynamics. An approach that has proven to be extremely successful is the Boltzmann-Uehling- 

■ Uhlenbeck (BUU) model of heavy ion collisions |4|. In BUU simulations, nucleons can suffer hard collisions and can 
CO . also move on curved trajectories owing to interactions with the self-consistent mean field. The interaction term we 

' use in this work is introduced and described in Refs. [^j6|. In connection with observables, a considerable amount of 
, information on the nuclear equation of state and on nonequilibrium precursor phenomena is accessible through the 

■ measurement of global collective behaviour. A large body of work in intermediate energy nuclear collisions has been 
' devoted to the measurement and to the theoretical calculation of both directed flow , and elliptic flow [^l]-^ . 

--^ . For reviews, see [Q, and references therein. 

1^ ■ Near the low end of the intermediate energy spectrum, some studies have put forward the possibility of observing 
I . experimental signatures of new phenomena. A good example is that of modified in-medium nucleon-nucleon cross sec- 
^ ' tions [p"5|-p0|. Those would signal a departure from vacuum properties. While the confirmation of such manifestations 
^ , would indeed be extremely interesting, one must keep in mind that such "new physics" issues must be addressed with 
Ph ' an approach that incorporates all of the known physics in a computationally tractable model. Our goal in this paper 
^1 is to first briefly describe such a model, and then to compare its results with experimental data. At the lower beam 
• ' energies considered in this work, the problem of energy and momentum conservation in transport models is a pressing 
/\ ' one. In addition, the momentum dependence of the nuclear mean field is an unavoidable feature, both from the point 
^ of view of theory and from that of experiment |2^-|2j]. With those issues in mind, we use a momentum-dependent 
- - ' lattice Hamiltonian model as a solution of the BUU transport equation. This specific model has previously been used 
in \2l], and provides for excellent energy and momentum conservation. A similar approach was developed in parallel 
in 1 26 ] . The momentum- independent lattice Hamiltonian method was first used in | |2^ ] . 

Our paper is organized as follows: the next section is a brief discussion of the momentum-dependent lattice Hamilto- 
nian model. The following section introduces the self-consistent in-medium correction to the free-space nucleon-nucleon 
cross section as dictated by the functional momentum dependence that is observed in the nuclear mean field. Next, 
a comparison of model results with compiled and recent measurements of directed and elliptic flow in the range of 
laboratory bombarding energy E^/A :25^800 MeV is performed. We then summarize and conclude. 



II. THE MODEL 



A. Lattice Hamiltonian solution of the mean field 



In this work, we use the BUU transport model which characterizes the time evolution of a system of nucleons during 
the course of a collision of two heavy ions. The BUU equation reads: 
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^IhhH +^.h- Vrfir, P, t)^Wph- ^pf{r, p, t) 
= j d^p, dn ' (^v;,i X ^) {f'fjh ffj'fi) . (2.1) 

The left-hand side of ( |2.l|) accounts for the time evolution of the semi-classical nucleon phase space distribution 
function / due to nucleon transport and mean field effects, and the right-hand side accounts for binary nucleon- 
nucleon coUisions which also modify /. The nucleon single-particle Hamiltonian h ~ t + u, contains both the single 
particle kinetic energy and the mean field single particle potential. The starred quantities in the collision term represent 
in-medium values of the relative velocity and differential nucleon-nucleon scattering cross section respectively, and / 
accounts for Pauli blocking Q . 

For the nucleon-single particle potential, we adopt the momentum-dependent parameterization for nuclear matter 
used in Q, supplemented with symmetry and Coulomb terms 
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This functional is amenable to computational applications and also incorporates nonequilibrium features 1^,^. The 
strong interaction part of the mean field is fixed by specifying the value of the parameters A, B, a, C and A, neglecting 
the isospin and Coulomb terms. The parameter set from ^J, deduced for symmetric nuclear matter, specifies: A — 
-322.0 MeV, B = 352.5 MeV, a = 12/11, C = 62.75 MeV and A = 1.58p'^ which gives u{p = 0) = -72.4 MeV, 
u(j) = pf) — —51.4 MeV and u{p oo) = -1-30.5 MeV for a zero temperature equilibrium nuclear matter distribution 
at saturation density po — 0.16 fm^"^. The zero temperature effective mass at the Fermi surface and nuclear matter 
compressibility for this parameterization are m* /m — 0.67 and K — 210 MeV respectively. The charge-dependent 
part of the nucleon mean field is taken into account with the isospin and Coulomb contributions to the single particle 
potential, where T3 is the third component of isospin of the nucleon and p„ and pp are the local neutron and proton 
densities, respectively. For the isospin part of the mean field, empirically it is found that D = 34±4 MeV. Furthermore, 
support for D ranging from 27 — 40 MeV is found from various phenomenological investigations |2^, and Refs. therein]. 
We adopt the value of D = 32 MeV previously used in the context of simulations of heavy ion collisions |29|-|3^ . The 
above parameterization of the mean field provides a good fit to the nuclear equation of state for both symmetric nuclear 
matter and pure neutron matter as calculated in p3| , p4| for a only slightly larger value of the nuclear compressibility. 
Note however that the exact high density behavior of the symmetry energy is still an open question that stems 
partly from the uncertainty in the high density many-body calculations [^, and partly from the different theoretical 
paradigms leading to the equation of state psj . 

The nucleon dynamics are realized through the use of the lattice Hamiltonian |2^ equations of motion, whereby 
the semi-classical nucleon phase space density for N^ns systems of A nucleons is projected onto a configuration space 
lattice with a finite grid constant 5x. At lattice site a, the semi-classical phase space density reads: 



/a(p,<) = ^ E R{rc.-m)Sip-m) (2.3) 



With this, the lattice Hamiltonian equations of motion for particle i read 
dn Pi{t) 



dt {pi{ty+m^y/ 



N,^,{8xf R(r^ - n{t))Vj;^u^{p„ T^t) (2.4) 



^ = ~N,ns{5xy2^U^{.P^,Tlt)V?^R{fo. - f,(i)). 



(2.5) 



where Ua{pi-,T^,t) is the discretized form of (2.2), and R is the finite width configuration space nucleon form factor. 



In this work, we use N^ns = 100, Sx = 1.00 fm and a quadratic form factor which provides for good energy and linear 
momentum conservation as shown in [ p6[ . In one dimension, the form factor is defined as R{x) — (3(5x^/4 — x'^)/6x'^ 
for < x < Sx/2, Rjx)^ {3Sx/2 - x)y2Sx^ for dx/2 <x < 36x/2 and R{x) = for 3Sx/2 < x. The normalization 
factor of 1/Nens in (2.3) and l/6x^ and l/2Sx^ in R are chosen such that particle number is conserved. That is: 
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at all times t. Note that the lattice Hamiltonian equations of motion (2.4) and (2.5) are derived from the total 
Hamiltonian iJ, and not the single particle Hamiltonian h^_as is the case with the often-used test-particle method [^p7[. 
This work constit utes the first application (with Refs. ||2^ and Pq]) of the lattice Hamiltonian approach with the 
interaction of Eq. (2.2) calculated fully dynamically. See also Ref. ]37|]. 



B. In- medium effects in nucleon collisions 



For the calculation of the Uehling-Uhlenbeck collision term, we follow the cascade procedure outlined in [Q in 
conjunction with a Pauli blocking algorithm which distinguishes between protons and neutrons. However, the in- 
medium scattering cross section is calculated such that the functional dependence on momentum that is observed in 
the nuclear mean field is respected |^^. Fermi's Golden Rule gives an expression for the total free space nucleon- 
nucleon cross section in vacuum: 

^rel 

where vf^i and o'^p are the relative velocity of the colliding nucleons and density of final states for the two-nucleon 
scattering process respectively. Both of these quantities are expressed as the free space value. In matter, we write the 
total in-medium nucleon-nucleon cross as 

where the starred quantities represent in-medium values. Detailed many-body calculations offer support for t*j^ tfi: 
as expected for low local densities and/or low energy collisions and Refs. therein]. We follow this line of thought 
here. In [^8| , the in- medium cross-section is calculated as it is here, however a simplified momentum-dependence in 
the mean field potential was used. We rewrite the in-medium elastic scattering cross section as: 

~* ^ ^'^ O-free- (2.9) 
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FIG. 1. Nucleon effective mass in equilibrium nuclear matter as a function of momentum at saturation density po (left panel) 
and density at the Fermi surface (right panel) for temperatures of 0, 10, 20 and 30 MeV. In the above, k — p/h. 



We note that a closed form solution to (2.9) exists for an equilibrium nuclear matter distribution. In this case, the 
in- medium nucleon-nucleon total cross section reads: 



m*{p,p) 



(2.10) 
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where the nucleon effective mass to* is defined as: 

'm*{p,p) ^ A ^ m du{p,p) 
m \ p dp 

and we have replaced the f dependence in u with a p dependence which is appropriate for an equiUbrium nuclear 
matter distribution. Note that the density p{f) is a constant in this case and it enters into the above relation self- 
consistently via the semi-classical phase space distribution function / (see, for example, Ref. p9{). In this work a 
modern and accurate parameterization of the nucleon- nucleon cross sections is used [ ^0[ . We show in figure |l| the ratio 
of the nucleon effective mass to the free mass as a function of momentum, density and temperature for an equilibrium 
nuclear matter distribution. As seen in the figure, the nucleon effective mass approaches the free mass in the high 
momentum, high temperature and low density limit as is expected from the many-body nature of this phenomenon. 



(2.11) 




FIG. 2. Ratio of the in-medium to that of the free space nucleon-nucleon elastic scattering cross section for an equilibrium 
nuclear matter distribution at T = MeV. 

In figure ^ we show the ratio of the in-medium to that of the free space nucleon-nucleon total cross section as 
a function of both density and momentum for an equilibrium nuclear matter distribution at T = 0. Note that for 
momentum-independent mean field potentials, the effective mass is equal to the free mass and consequently, the 
in-medium cross section is equal to the free space cross section. 



C. Dynamical effects in heavy ion collisions 



We now wish to apply our model developed thus far to simulate the collision of two heavy ions. In this scenario, 
dynamical effects will generate a nuclear matter distribution that is in general not that of a zero temperature system. 
Furthermore, one also expects a nuclear matter distribution that is different fro m th e equilibrium situation. In th is 
case, one cannot use the zero temperature equilibrium closed form solution of (2.2) as given in g. Instead, (U), 
(2.4), (2.5) and ( |2.9D must be explicitly and dynamically calculated at all times during the collision process. In this 
work, we present the first results of the self-consistent calculation of ( |2.9| ) as applied to the dynamical collision of 
heavy ions. 

In the initial stages of the collision of two heavy ions, large deviations from equilibrium will result in binary nucleon- 
nucleon centre of mass collision energies higher than that at later stages in t he co llision process. From figure || one 
can infer that initially, owing to the density and momentum-dependence in (2.1C), the in- medium cross section will 
be larger at earlier stages in the collision compared to later times. In addition, one can also infer that higher centre of 
mass nucleus-nucleus collision energies will result in an in- medium cross section that is closer to the free space value. 
We illustrate this for a simulation of the head-on collision of two ^"^Bi nuclei at laboratory bombarding energies of 
25, 150, 500 and 1000 MeV/A in the left panel of figure |. In this figure, we show the mean value of the in-medium 
nucleon-nucleon cross section for unblocked collisions as a function of time. We note that previous attempts to model 
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the in-medium nucleon-nucleon cross section have used a constant overall scaling factor |15yiC 
a phenomenological density-dependent in-medium cross section of the form 



Others have used 
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FIG. 3. Mean value of the ratio of the in-medium to that of the free space nucleon-nucleon scattering cross section for 
dynamical collisions of '^"^Bi-l-^^^Bi at laboratory bombarding energies per nucleon of 25, 150, 500 and 1000 MeV as a function 
of time. The left panel displays the results of the self-cor isiste nt calculation of ( ^.jj] ) and the right panel shows the results of 
the phenomenological density-dependent ansatz given in ( ^.12| ) with q = 0.20. 

cr* = (1 - ap/ po)(Jfree, (2-12) 

where a is an adjustable parameter We show in the right panel of figure || that adoption of this simple 

ansatz results in an in-medium cross section that has considerably different behaviour from that of the self-consistent 
calculation. The self-consistent modification to the free space cross section in ( |2.9D effectively probes the shape of the 
momentum-dependent part of the single particle potential at all momenta. Such a detailed modeling is not possible 
with neither a constant scaling factor nor the density-dependent relation in ( 2.12 ). We will address the consequences 
of these differences in terms of elliptic flow in section II] . Before we leave this section, we show in figure ^ a sample 
(250 events at each energy) of the values of the self-consistent in-medium cross section as a function of density for the 
t : 10 ^ 20 fm/c time slice from figure ^ for all energies. 
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FIG. 4. Value of the self-consistent in-medium cross section for the t : 10 ^ 20 fm/c time slice from figure ^ for all energies 
considered in that figure (scattered points). Moving from left to right and top to bottom, the panel s are for energies of 25, 150, 
500 and 1000 MeV. The solid lines represent the value of the in- medium cross section using (2.12) with various values of the 
parameter a. 
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Owing to the additional momentum-dependence in the in-medium cross section, we see that for a given density, 
the former can take on drastically different values. The solid lines i n tha t figure indicate the value of the in-medium 
cross secti on that would be obtained by using the relation given in ( ^.12 ) for two different values of a. Again, we see 
that ( 2.12 ) clearly oversimplifies the in-medium result from (2.9). 

To summarize, we note that the calculation of the self-consistent momentum (and density) dependent in-mcdium 
cross section indicates that the latter approaches the free space cross section for high centre of mass nucleus-nucleus 
collision energies and steadily decreases with centre of mass nucleus-nucleus collision energies. Also, at early stages 
in the nucleus-nucleus collision, non-equilibrium effects results in an in-medium cross section that is closer to the 
free space value then compared to later times in the nucleus-nucleus collision process. This implementation of a 
self-consistent, momentum-dependent, BUU transport model represents a considerable numerical challenge. Efforts 
to develop algorithms in a parallel architecture are under way. 



III. COLLECTIVE EFFECTS 



A. Directed Flow 



One of the often discussed signatures of collective motion in heavy ion collisions is the directed flow [Q . Typically, 
the projection into the reaction plane of the transverse momenta of the reaction products is plotted against the 
rapidity, resulting in the well known "S-shaped" plot. The directed flow is obtained by taking the slope of this plot 
at mid-rapidity. Specifically, one fits the flow plot with: 



d <Px > 
dy 



ai + a2y + a^y 



(3.1) 



where the rapidity fit interval is tak en t o be \y\ < yproj, and yproj is the initial projectile rapi dity. In this case, the 



flow F is given by the value of 02 in (3J). Often, a narrow mid-rapidity region is selected, and ( p.l[ ) is used to fit the 
flow, with 03 0. Some recent flow data from the FOPI collaboration is shown in figure ||. 
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FIG. 5. Directed flow F^°'^ as a function of impact parameter for Au+Au collisions at 400 MeV/A. The measured points |42[| 
shown by triangles, are obtained using a coalescence model as described in the text. T he r esults from the BUU calculations are 
shown by the solid and open circles for a free space and in- medium cross section from (2^) respectively. The BUU calculations 
are shown for all protons at mid-rapidity. The curves are drawn to guide the eye. 



In this figure, the directed flow for Au-|-Au collisions at 400 MeV/A is shown as a function of impact parameter. The 
experimental data in this plot was obtained by using a coalescence model which weights all detected fragments by its 
measured charge . We also show in this figure the re sults of our BUU simulations with a free space nucleon-nucleon 
cross section, and the in-medium cross section given in (2.9). Here, the normalized flow = F x {yproj /Pproj) cm is 
used. The coalescence-invariant prescription was develo ped to enable direct comparisons with transport models such 
as the BUU. All the BUU results were obtained by using (3T) with and without 03 = 0. For the = \y/yproj\ < 0.45 
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cut, we find the two methods to be within one standard deviation of each other. Qualitatively, we find that both the 
free space and in-medium simulations predict the observed maximum in the flow magnitude for intermediate impact 
parameters. More quantitatively, we find the free space cross section over-predicts the flow magnitude for all impact 
parameters. The calculations with the in-medium cross section however, are generally in good agreement with the 
data at all impact parameters. 
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FIG. 6. Directed flow excitation function for Au-|-Au collisions. The measured points shown by the triangles and stars 

are for fragments of charge Z = 1,2. The measured points shown by the diamonds are for fragments of charge Z — 2. The 
BUU results (circles) are as in figure ^ and are impact parameter averaged from 3.5 fm< b < 6.5 fm for Ek/A > 100 MeV at 
rapidity —0.19 < yn < -1-0.29. For E^/A < 100 MeV, b < 0.39 at rapidity |j/„| < 0.5, fm as is also the case with the measured 
points shown by the diamonds. The curves are drawn to guide the eye. 

The directed flow excitation function also provides a means by which comparisons of measured signals to model 
calculations may be performed. With this in mind, we show in figure |^ a compilation of data near the balance 
energy for central impact parameters and charge z = 2 fragments and for semi-central impact parameters at 
Ek/A : 150 — > 800 MeV for charge z = 1,2 fragments. The BUU results for the central impact parameters 

are shown for |?/„| < 0.5 and b — h/b^ax < 0.39 fm and for —0.19 < y„ < -1-0.29 and (3.5 < b < 6.5) fm for the 
semi-central impact parameters. To cornpare our BUU results with the measured fragment flow, free nucleons are 
omitted from the flow analysis; see Ref. for a discussion. 
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FIG. 7. Same as figure y, but for free protons only. The measured points ||9| are shown by triangles and the BUU results 
(circles) are as in figure n. The curves are drawn to guide the eye. 
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For the entire energy range probed here, we find that the calculations employing the in-medium nucleon-nucleon 
cross section result in a reduced value of the directed flow as compared with the free space cross section. From this 
figure, we find the in-medium result to be in better agreement with the data. A more stringent test of the model in 
terms of the directed flow is shown in figure 0. In this figure, the data for free protons ||] is shown along with the 
model predictions for free protons. The figure indicates that the predicted directed flow over-estimates the observed 
value. We point out however, that the coalescence-invariant comparison in figure ^ with the more recent data is in 
better agreement with the model calculations. 



B. Elliptic Flow 



We now turn to comparisons of the measured elliptic flow with our BUU model. The elliptic flow is a measure that 
quantifies the azimuthal anisotropy of the momentum distribution. Specifically, we fit the azimuthal distribution of 
nucleons about the reaction plane with a Fourier expansion of the form: 



dN 

TJT = Po (1 



2picos{(f>) + 2p2Cos(2(f))) 



(3.2) 



where the ellipticity coefficient p2 depends on the in-plane and out-of-plane flow amplitudes. In addition, the 
anisotropy ratio is deflned by: i?„ = (1 — 2p2)/{l + 2p2). A ratio i?„ larger than imity signals a preferred out- 
of-plane emission. 
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FIG. 8. Mid-rapidity (|j/n| < 0.15) free-proton anisotropy ratio for Bi-|-Bi collisions at incident laboratory bombarding 
energies of E^/A =200, 400, 700 and 1000 MeV moving left to right and top to bottom. The measured points are shown 
by the squar es. T he BUU results are for a free space cross section (solid circles), den sity-dependent in-medium cross section 
(stars) from (2.12) and self-consistent in-medium cross section (open circles) from (2.9). The impact parameter for all cases is 
b ~ 0.65. There are no data points available at 200 MeV. The curves are drawn to guide the eye. 



In it was shown that BUU calculations of the elliptic flow for a Bi-|-Bi system with a momentum-dependent 
nuclear mean field favoured a nucleon effective mass of r a* /m ~ 0.6 — 0.7. We show in figure ^, the results that we 
obtain for BUU calculations with the mean field in ( ^.2D of effective mass m* / m ^ 0.67, for b oth a free space cross 
section and separately, the in-medium cross sections obtained with equations (2.9) and (2.12). The figure indicates 
that adoption of the self-consistent correction to the free space cross section is significant for the anisotropy ratio at 
the lower energies only (200 and 400 MeV/A). For the higher bombarding energies, the self-consistent result does not 
differ significantly from the free space result. Since the in-medium effect from (^^) probes the shape of the momentum- 
dependence in the mean field it is more effective for lower bombarding energies, with many nucleon-nucleon collisions 
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occurring near the Fermi surface where the momentum-dependence is steepest. At high bombarding energies, the 
nucleon-nucleon coUisions occur far from the Fermi surface where the momentum-dependence in the mean field is 
relatively flat, and as such, the correction to the free space cross section is small. The experimental results are from 
p6| . Unfortunately there are no data at the lower energy, where the differences in theoretical results are the largest. 
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FIG. 9. Excitation function for the ellipticity coefficient for Au-|-Au collisions. The measured points shown by the 
triangles are for fragments of mass A = 1,2, 3, and are for an impact parameter range of b : 5.3 ^ 7.3 fm. The BUU results 
(circles) are as in figure |^ and are shown for an impact parameter of b = 6.5 fm for all mid-rapidity nucleons. Mid-rapidity 
particles are selected with a < 0.1 cut. The curves are drawn to guide the eye. 
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We now turn to the elliptic flow excitation function for Au-I- Au recently measured by the FOPI collaboration . 
In figure ^, we show the ellipticity coefficient measured for fragments with A = 1,2,3 along with the appropriate 
BUU result. W e find that both parameterizations of the nucleon-nucleon cross section (free space and in- medium 
result from (^^)) are able to qualitatively reproduce the experimental trend. Quantitatively, the calculation with the 
in-medium cross-section is better able to reproduce the measured signal. 
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FIG. 10. Same as figure |9| but for high transverse momentum protons only {p^ > 0.8). The measured points are taken from 
and are for an impact parameter range of b : 5.3 7.3 fm. The BUU results (circles) as in in figure |^ and are shown for 
free protons at b = 6.5 fm. Mid-rapidity particles are selected as in figure H. The curves are drawn to guide the eye. 



We note that the ellipticity parameter is calculated in a frame where the principle axis of the momentum ellipsoid 
coincides with the beam axis. A more stringent test of the model is shown in figure fol where the ellipticity coefficient 



9 



for free protons with transverse momenta p'^' >0.8 only, is plotted as a function of incident bombarding energy. Note 

that p"^ = p^/pbeam- From this figure, we find that the elUptic flow in the energy range (150 < Ek/A < 600) MeV 
is over-predicted by the model. This agrees with figure |[ where we see that the magnitude of the elliptic flow is 
over-predicted at high transverse momentum. This identifies the phase space areas where the BUU seems to perform 
perhaps less satisfactorily. However there, the phase space distribution functions are small and therefore those regions 
do not contribute importantly to the momentum-integrated observables. This is seen in our results as well as in the 
experimental data. This observation however represents valuable information for future theoretical developments. 

IV. CONCLUSION 

In summary, we have shown for the first time the effect of the parameter-free self-consistent calculation of the 
nucleon-nuclcon in-medium scattering cross section in BUU calculations of both directed and elliptic flow from 
Ek/A : 25^800 MeV. Our results indicate that, globally, the in-medium cross section better describes the flow data 
than with the free space cross section. In addition, the self-consistent calculation of the in-medium cross section results 
in a value of the latter that is less than the free space cross section for low bombarding energies and approaches the 
free space cross section at high bombarding energies. While in some sense transport models utilizing vacuum cross 
sections are a well-defined low-density approach, there is still some work to be done to formulate and use a totally 
self-consistent many-body transport theory. For example, one needs to generalize to the coupled set of transport 
equations that include all additional degrees of freedom |48] . The work in this direction is continuing. 
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